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Abstract 

We have derived the popularly used parametrization formulae for quark masses at 
low densities and modihed them at high densities within the mass-density-dependent 
model. The results are applied to investigate the lowest density for the possible 
existence of strange quark matter at zero temperature. 
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The possible existence of strange quark matter (SQM) has been a focus 
of investigations [1] since Witten’s conjecture [2] that quark matter with 
strangeness per baryon of order unity may be bound. Because of the well- 
known difficulty of Quantum Chromodynamics in the nonperturbative do¬ 
main, phenomenological models reflecting quark conhnement are widely used 
in the study of hadron, and many of them have been successfully applied to 
investigate the stability and properties of SQM. One of the most famous mod¬ 
els is the MIT bag model [3] with which Farhi and Jaffe [4] hnd that SQM 
is absolutely stable around the normal nuclear density for a wide range of 
parameters. 

Another popularly used model is the mass-density-dependent model in which 
quark conhnement is achieved by requiring [5] 

lim m = oo, (1) 

rib^O 


where m is the quark mass, Ub is the baryon number density. By using this 
model, Chakrabarty et al. [6] obtain a signihcantly different result: only at 
very high densities does SQM have the absolute stability. 

Benvenuto and Lugones [7] point out that this is the consequence of an in¬ 
correct thermodynamical treatment. They add an extra term to the energy 
expression, and get similar results to those in the bag model. 

S. Chakrabarty et al. [6,8] have already discussed the limitation of the con¬ 
ventional MIT bag model which assumes that the quarks are asymptotically 
free within the bag. In order to incorporate the strong interactions between 
quarks, one has to fall back on the perturbative theory, while the mass-density- 
dependent model mimics not only the quark conhnement, but also the inter¬ 
actions between quarks. Therefore, the validity of quark mass dependence on 
density is of utter importance. 

However, the popularly used quark mass formulae [5-7] are pure parameter- 
izations without any real support from underlying held theories up to now. 
Our motivation in writing this letter is to try to derive the relation between 
the quark mass and density from more fundamental principles on one hand, 
and to study the critical density for SQM on the other hand. By “critical 
density” we mean in this paper that only above the density does SQM have 
the possibility of existence. 

In a recent work [9], we have demonstrated that the quark mass and quark 
condensates satisfy the following relation 

m 1 

me 1 - (gg)n6/(gg)o’ 


( 2 ) 
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where we have suppressed the color index, the model parameter is changed 
to nic from nigo which is saved for the original mass, {qq)^ and {qq)nt, are 
the quark condensates, respectively, in vacuum and in strange quark matter 
with baryon number density rib. In this paper, we will see that Eq. (2) has 
important properties of both conhnement and asymptotic freedom. 


According to the obvious equality 


lim 


{m)n, 

{qq)o 


1 , 


( 3 ) 


we can expand the relative condensate as 

(W)n, T . , 

= 1-h higher orders m n;, + ■ • •, 

m)o « 


( 4 ) 


where 


/ d {qq)n, \ ^ 
\dnb {qq)o / 


( 5 ) 


At low densities, we can naturally ignore all terms in Eq. (4) with orders in 
rib higher than 1 and obtain 


{QQ)ni, ^ ^ 

{qq)o a'' 


Substituting Eq. (6) into Eq. (2), we get m = rrica'/rib. 

At high densities, the approximation (6) is no longer valid. We will soon see 
that the quark mass is inversely proportional to n^/ , rather than just rib. 

The Hellmann-Feynman theorem [10] indicates that 

(*(A)| A//(a)|4.(A)) = T(4,(a)|//(A)|>I-(A)), (7) 

where H{X) is a Hermitian operator, |'h(A)) is the normalized eigenvector of 
H{X), X is an independent real parameter. 

We express the effective Hamiltonian density as 

hfeff = H' + mqq, (8) 
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where H' is the kinetic energy term. In the quark mass-density-dependent 
model, the total (or effective) mass depends on the density. It can be divided 
into two parts, namely, m = tuq + mj, where rriQ is the original mass indepen¬ 
dent of density, and mi is the interacting mass mimicking the strong interac¬ 
tion between quarks. Assuming that the quark condensate depends merely on 
the total quark number density, or, is blind to the conhguration of the system, 
we, for simplicity, consider only one flavour case in the following derivation. 

Substituting / (PxH^q for H{\), and mo for A in Eq. (7), we have 

(9) 

Because the expectation value of ifgfj should equal that of T^qcd, we write 
down 


Applying this equation respectively to the cases I'h) = \nb) and |\k) = |0), 
where \nb) denotes the ground state of the quark matter at rest with baryon 
number density Ub, |0) is the vacuum, and then taking the difference, we obtain 

(gg)n, ^_ ]__d^ _ 

{Q(1)o \qq\odmo dmo \qq\odm' 


where the uniformity of the system has been taken into account, \qq\o = 
— {qq)o, e is the energy density of the quark matter. Compared with the cor¬ 
responding formula for nuclear matter [11], they are very similar in form. But 
the physical contents are very different. The energy density there depends 
on the nuclear density and the quark current mass which is an independent 
quantity. Here we have the following chain relation 


/nb 

^mo 

Amd .mo 

^mj/ 

^Ub- 


This is why we should take mo, rather than m, as the substitute for A. 
Combining Eqs. (2) and (11), we obtain 

de 

m—= mc\qq\o = B. (12) 
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At zero temperature, one has already known 


e 


4 

gm 


x{2x^ + 1) + 1 


sh , 


(13) 


where sh = ln(a; + vx^ + l), x = Pf/m, the Fermi momenta pf = 
g is the degeneracy factor(6 for quarks). If replacing the e in 
Eq. (11) by this expression, we can see that the relative condensate indeed 
has the similar expansion to Eq. (4). 

Substituting Eq. (13) into Eq. (12), we have 

xV+ 1 — sh~Vx) = ^ ■ (14) 

gm^ 


Generally, this equation can be solved numerically. However, It can be ana¬ 
lytically solved both at low and high densities. 


At low densities, we expand the left hand side of Eq. (14) to x^ term, and get 



1 

2 


1 / ^ 31 ^ y/^ 

2 ^ 5^2/352 J 


(15) 


If we only expand Eq. (14) to x^ term, or equivalently, ignore the second term 
in the parentheses of Eq. (15), then 


B 



(16) 


This result was obtained from other arguments by the previous authors [5] 
many years ago. It agrees to the conhnement property of quarks, and has 
been popularly applied to investigate the properties of SQM [6-8]. 


At high densities, Eq. (14) becomes 




47r^i? 

gm? 


(17) 


The hrst term on the left is much greater than the second one. So we ignore 
the second term and get 


m = 


2n^Bjg 

Pf 


(47r/9/V^)V3v/]e 



(18) 


5 



A more precise solution can be obtained by substituting this expression for 
the m's in the second term on the left hand side of Eq. (17): 


m = 


2TT^W/~g 


Pf 


l + ^ln(^ 


4 

Pf 


7r‘^B/g' 


(19) 


Equation (18) or (19) is in accordance with the asymptotic freedom of QCD. 
Table 1 

Comparison of the relative errors(%) for different mass formulae. The densities are 
in unit of the normal nuclear density uq = 0.17 MeV-fm“^. Parameter (3 is taken to 
be (137 MeV)^. means over 100%. 
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* 

Eq.(16) 

0.40 
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Eq.(18) 
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3.23 

1.53 

0.98 

Eq.(19) 
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43.09 

10.75 
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0.44 

0.23 

0.15 


In Table 1, we show the relative errors for different mass formulae. We see 
that the quark mass is inversely proportional to density at very low densities. 
When the density exceeds the nuclear density, the logarithmic term occurs, 
and the quark mass changes asymptotically to be in reverse relation to 
Therefore, we should apply Eq. (18) to the investigation of SQM, namely. 


^u,d — 


rus = 



p, 



( 20 ) 

( 21 ) 


The parameters (3q and (3s determine the critical density Uc rather simply 


Uc 


- Plf' 

ir^/7 


( 22 ) 


where 7 is dimensionless: 1.5 < 7 < 2. We use 7 = 1.995 in the calculations. 

To prove this expression, we assume the SQM to be a Fermi gas mixture of n, 
d, s quarks and electrons with chemical equilibrium maintained by the weak 
interactions: d, s u + e + W, -s + n <->• n + d. At zero temperature, the 
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thermodynamic potential density of the species i is 




dSTT^ 






m‘f +3mj In 






mj 


rrii 


. (23) 


where i = u,d, s, e(ignore the contribution from neutrinos); Qi is the degener¬ 
acy factor with values 6 and 2 respectively for quarks and for electrons; mu,d 
and rus to be replaced by Eqs. (20) and (21). 

The corresponding baryon number density is 

-It - 


For a given Ub, the chemical potentials /ij(i = u, d, s, e) are determined by the 
following equations [4] 


P'S — h) 
fJ'u T h'e h) 

Ub = -{nu + nd + Us), 

2 11 

-n„ - -rid - -Us - ne ^ 0. 


(25) 

(26) 

(27) 

(28) 


The last two equations are equivalent to 


nd + ns + ne = 2nb. 


(29) 

(30) 


We therefore dehne a function of fie 


F{fie) = (/i^ - + (/i^ - + -{fil - - 27rV, (31) 


where 


fi 


he + 



+ Wnb +-{fil 


m2) 3/2] 2/3. 


(32) 


Because mg > mu,d, the equation F{fie) = 0 for fie has solution if and only if 
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Fifie) < 0 . 


(33) 

(34) 


At the critical density Uc, the equality signs in the above should be taken, and 
so we can easily hnd 


l.STT^nc < {ml — < 27r^nc. 


(35) 


Solving for Uc from this inequality, we accordingly obtain Eq. (22). 
For a more exact Uc, we can numerically solve the following equation 


ms - \Jml^^ + [?>'K^nc - {ml - ,^) 3 / 2 ] 2/3 

= 'Jml + {3[27r2nc — {ml — 
in the range 

~/27, ‘ VlE^ 


( 36 ) 


( 37 ) 


In fact, the expression (22) is precise enough for practical applications because 
of the smallness of electron content. 

At Tic, the strangeness fraction becomes zero. When the density decreases 
further, the equation group (25-28) which determines the conhguration of the 
system has no solution. This indicates that Uc is the lowest density for the 
possible existence of SQM. 

To calculate the critical density, we need to know the concrete values of Pq 
and Ps- 

Generally, there are two constraints on Pq [4,6,7]: firstly, it should be so big 
that the energy per baryon for light-flavour quark matter is greater than 930 
MeV in order not to contradict standard nuclear physics; secondly, because 
we are interested in the possibility of absolute stability of SQM, it should not 
exceed an upper limit so that the energy per baryon can be less than 930 MeV 
for symmetric three flavour quark matter. Therefore, Po must be in the range 
from a hxed lower limit Pomin to an upper limit Pomax- By using the method 
in Ref. [7], we obtain: /3omm = (135 MeV)^ Pomax = (147 MeV)^. 

Because at very low densities the mass formula should be Eq. (16), we apply 
this formula to estimate the lowest permissible value for B independently by 
the same method. The result is that B should be no less than 69 MeV-fm“^. 
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This requires /9o to be greater than (137 MeV)^ according to the relation 
(3 = (47r/9/^)i/3/B (see Eq. (18)). It is remarkable that the requirement is 
in accordance with the lower limit of Pq. With our purpose of studying the 
lowest density for the possible existence of SQM in mind, we take Pq = (137 
MeV)2. 

As for j3s-, there is no reliable method to £x it presently. If regarding the 
nucleon as a special SQM with zero strangeness, one has 



where the ratio Us/<yN has been taken to be 6 [12]. The corresponding ric is 
2.7no, not too far away from the normal nuclear density. Therefore, wether 
SQM can exist near the nuclear density is still an open question to study. 
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